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Abstract. We consider the totally asymmetric simple exclusion process (TASEP) 
on the periodic chain in the presence of a single impurity site that is inaccessible to 
other particles and therefore acts as a static defect. Particles are allowed to advance 
any distance I > 1 on the right with the probability that decays as l~( 1+a \ where 
a > 1. Despite the long range of hopping, we find the same type of phase transition 
that occurs in the standard short-range TASEP with a defect site where defect induces 
a macroscopic shock in the stationary state. In particular, our model displays two 
main features characteristic of the short-range TASEP with defect site: a growth of 
the shock width with system size L as L 1 / 2 or L 1 / 3 , depending on the existence of the 
particle-hole symmetry, and the power-law decay in density profiles of the shock phase. 
However, unlike the profiles in the short-range case, we find that the latter are well 
reproduced by the mean-field approximation, which enables us to derive the analytical 
expression for cr-dependent exponent v = a — 1 of this power-law decay and the point 
er c = 4/3 at which the transition takes place. 



PACS numbers: 02.50.Ey,05.60.-k,47.40.Nm 
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1. Introduction 

Driven non-equilibrium systems have been a subject of extensive study in the past 
few decades. Among many others, systems of interest include growing interfaces pQ, 
vehicular traffic [2 J and biological transport [3]. One of the simplest of related models is 
the totally asymmetric simple exclusion process (TASEP) [I], in which particles jump 
to the nearest-neighbour site on the right, while each site holds at most one particle. 
In case of the open boundary conditions, the finding of the exact stationary state 
[51 [6] exhibiting boundary-induced phase transitions [7J has triggered a study of various 
generalizations including Langmuir kinetics [HI E] or inhomogeneous hopping rates [10J. 
The last example is of a particular interest in studying the effects of disorder (see, for 
example, [H]), which is hard to avoid in most realistic situations. 

Disorder in TASEP is usually introduced through inhomogeneous hopping rates 
that are either associated to particles or sites. Due to the long-range correlations already 
present in the pure TASEP, strong effects are found even in the presence of a singe defect 
particle |L2] or site [131 HH H51 EE EL71 HH] • In the case of a single defect particle, the exact 
solution was found using the matrix-product ansatz [12] and was recently generalized 
to arbitrary number of defects [19]. On the other hand, TASEP with defect sites (and 
even a single one) is still an open problem, usually dealt with Monte Carlo simulations 
and various mean- field approaches (see [20J and references therein). 

There is, however, a particular generalization of TASEP which reproduces its main 
features, but is remarkably well described by the mean-field approach. Originally 
motivated by the emergence of long-ranged power-law-type correlations in TASEP, a 
generalized model was proposed [2TJ in which particles are allowed to advance any 
distance I > 1 with the probability that decays as l~( 1+a \ In case of the open boundary 
conditions, albeit modified due to the long range of hopping, this generalized model 
has the same phase diagram as the short-range case consisting of the low-density, high- 
density and the maximum- current phase, but with different effects at the transition lines. 
Careful analysis of this model [22] reveals the emergence of an effective bulk reservoir 
(similar to the one introduced by Langmuir kinetics), which for 1 < o < 2 dominates 
over the fluctuations. 

Similarity with the short-range TASEP and the applicability of the mean-field 
approach motivates us to study the effect of a single impurity also in the long-range 
case. At first sight, one could expect that a localized defect in system with long-range 
hopping would not have any significant effect on the flow of particles. The main objective 
of this paper is to show that this is not the case. For that purpose we introduce the 
simplest defect in terms of a static impurity that occupies one site on the lattice and 
is not allowed to hold any particle. This induces the same type of transition as in 
the short-range case, but with respect to a as a control parameter. We show that the 
model is well described by the mean-field approach which enables us also to obtain the 
analytical expression for a c at which the transition takes place. 

The paper is organized as follows. In section [2] we define the exclusion process 
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with long-range hopping on the periodic lattice and introduce a defect site to study its 
effects on the stationary states. Density profiles and their scaling properties obtained by 
Monte Carlo simulations are displayed in section [37X1 and analyzed within the mean-field 
approach in section 13.21 A brief summary of results is given in section HI 

2. Exclusion process with long-range hopping and impurity 

In the pure long-range model without defect [2Tj . N = pL particles are distributed on 
L sites of a one-dimensional lattice with periodic boundary conditions such that each 
site is either occupied by a particle [r n = 1) or empty (r n = 0). During an infinitesimal 
interval dt, a randomly chosen particle at site 1 < n < L attempts to jump I sites to 
the right, where distance 1 < I < L is chosen according to the probability distribution 
pi = /~( 1+cr ) / Xl(o~ + 1) with being the partial sum of the Riemann zeta function. If 
the target site n + I is empty, the move is accepted; otherwise, it is rejected. Compared 
to the standard (short-range) TASEP, common features include factorized steady state 
and the particle current of the same form, j = \(o~)p(l — p), where the additional factor 
^( <J ) = J2i I ' Pi = Cl(c) / '(l(o~ + 1) m &y be interpreted as the average hopping length. 
The model is well defined for a > 1 where A (a) is finite. For 1 < o < 2 it displays true 
long-range behaviour, while the short-range regime sets in for a > 2 [211 122] . 

To study the simplest case in which the flow of particles is obstructed by a single 
defect, an additional site is included at the position L + 1 and the periodic boundary 
conditions are shifted accordingly (r„ = r n+L+ i). The additional site is not allowed to 
hold any particle, but the particles are still able to pass it as long as a < oo (the flow 
is fully blocked only in the short-range limit a — > oo where p\ = <^,i). The key benefit 
from such definition of impurity is that it obstructs the flow of particles but does not 
introduce any new parameter in the model. Any change in behaviour therefore must 
come solely from the variation of a. 

3. Density profiles 

3.1. Monte Carlo simulations 

Typical density profiles obtained by Monte Carlo simulations for p = 1/2 are shown in 
figure CD As in the short-range case, one observes two regimes (separated by the value a c 
of the control parameter): for a > a c the defect is strong enough to induce shocks with 
densities p + and p_ = 1 — p + left and right of the impurity, respectively; for 1 < a < o c 
the density profile is shock-free. In the latter case the profile organizes itself in a way 
that maximizes the current, which attains the maximal value as in the pure long-range 
model, A(er)/4. 

For any finite system the shock is also characterized by its width, as a result of 
the fluctuations in the microscopic position of the instantaneous shock that performs 
random walk around the lattice. Using the useful concept of the second-class particle 
[T3] it was found that the width in the short-range model scales with L as L 1//2 and L 1 / 3 
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Figure 1. Density profiles in the (a) shock- free phase for a = 1.1 and (b) shock phase 
for a = 1.8 obtained by Monte Carlo simulations for L = 6400, p =1/2 and i = 10 8 
MCS/site. 



depending on whether p 7^ 1/2 or p = 1/2, respectively, which is related to whether 
the collective velocity v g = 1 — 2p at which the density fluctuations travel is finite or 
vanishes [I 3|. 

We recover this result also in the long-range model, in accordance with our earlier 
result [22] that the density fluctuations in the infinite system travel with the velocity 
v g = A(cr)(l — 2p). In figure [2] we plot the width £r,, obtained by the Gaussian fit of the 
discrete derivative (r n+1 ) — (r n ) vs. system size L for densities p = 0.5 (figure^) and 
0.55 (figure Eb). 
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Figure 2. Shock width £l obtained from the stationary density profile for various a 
and densities (a) p = 0.5 and (b) p = 0.55. Dashed lines have slopes (a) 1/3 and (b) 
1/2. 



Another interesting feature to compare with the short-range model are the site- 
dependent corrections to the bulk densities p_ and p + . In the short-range model it was 
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found that they exhibit power-law behaviour away from the defect site 

(r„) - p- ~ -n~" ', (tl-u) - P+ ~ 1 < n < L, (1) 

while different values of v were given in literature. Early results of Janowsky and 
Lebowitz for the periodic boundary conditions proposed that v — 1 [13] based on the 
Monte Carlo simulations. The problem was later studied in the context of surface growth 
[T6] in the particular model that can be mapped directly to TASEP with one defect site. 
They considered only the shock-free phase and obtained the exponent v = 0.41 using 
renormalization-group-like treatment of the exact solution for N = 2 and the Monte 
Carlo simulations. These results were later revised by Ha et al. [18] based on numerical 
results and analytical arguments arguing that the power-law behaviour ([1]) occurs in 
both phases, with v — 1/2 in the shock phase (p_ < p + ) and v = 1/3 in the shock-free 
phase (p_ = p + = 1/2) where the current is maximal. 

In order to check the power-law behaviour (TjQ) in our model, we assume the following 
scaling relation for a deviation of density from its bulk value (tlu) 

Ap_(n,L) = (r n ) - (r L/4 ) = L^f^n/L). (2) 
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Figure 3. Deviation of a density profile from its bulk value (tl/ 4 ) in the low-density 
domain of the shock phase for (a) a = 1.6 and 1.8 and (b) a = 2.0 and 2.3. Profiles 
have been scaled for different system sizes 3200, 6400 and 12800 according to ([2|) with 
(a) v = 0.6 and 0.8 and (b) v = 1, respectively (p = 1/2 and t = 10 8 MCS/site). 

Scaled density profiles are shown in figure[3h. for o = 1.6 and 1.8 and various system sizes. 
The best collapsing fit is achieved for v w 0.6 and v w 0.8 respectively, suggesting that 
v is equal to a — 1. On the other hand, for o > 2 one generally expects the short-range 
regime to set in [22] with a a-independent v. Indeed, figure [3b for a = 2 and o = 2.3 
shows good agreement between the scaled profiles if one assumes that v(o~) = u(2) = 1. 
This exponent matches the one suggested by Lebowitz and Janowsky in the short-range 
TASEP [13J. Let us also note here that in the shock-free phase the scaling is not so 
clear and remains elusive. Although we observe density profiles that certainly exhibit 
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long-range correlations, they are not described by the single exponent v as in We 
shall return to this point later in the discussion of the mean-field approach. 

Finally, let us touch upon the value of a c . The main problem in finding o~ c 
numerically are the finite-size effects that may create shock for any o provided that 
the system size is small enough. Near o~ c , system sizes needed to observe the transition 
become inaccessible to computer simulations and one has to relay on other methods. For 
example, in the short-range case where the transition takes place at value r c , Janowsky 
and Lebowitz suggested the estimate r c > 0.8 based on the exact expression for the 
current obtained for small systems [23]. On the other hand, Ha et al. [H] assumed 
the finite-size scaling of the order parameter p+(r, L) — p~(r, L) and obtained the best 
collapse of data for r c = 0.80(1). Here we take a different and less precise approach by 
looking at the smallest a for which the slope of the shock in the middle becomes steeper 
as one increases the system size L. In other words, for a < o~ c increasing L will flatten 
up the shock, while for a > a c one expects convergence to the step-like function. For 
the largest system size that we used in Monte Carlo simulations, L = 12800, this gives 
the estimate 1.32 < a c < 1.4 (figure Hj). 
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Figure 4. Emergence of the shock-free phase (a) for a = 1.32 and shock phase (b) for 
a = 1.4 from the density profiles on large systems, L = 6400 and 12800 (p = 1/2 and 
t = 10 s MCS/site. 



3.2. Mean- field approach 

The mean-field approximation, in which one neglects correlations by replacing (r n r m ) 
with (r n )(r m ), has been the starting point for many of the analytical approaches to the 
problem of defect site in TASEP (with random- sequential update) so far. The reason 
for this lies in the inapplicability of the usual tools like the matrix-product ansatz that 
proved to be extremely useful in the pure TASEP. Here we take the same approach of 
neglecting correlations between sites, which unlike the short-range TASEP gave accurate 
results in our earlier study of the pure long-range case [22] . 
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Before writing down the lattice equations for density profile (r n ) and applying the 
usual mean-field approximation (r n r m ) — > (r n )(r m ), let us recall the simple mean-field 
(SMF) approximation dealing only with bulk values p3] to estimate p_ and p+. In 
particular, we neglect any correlations between sites and assume the following density 
profile in the shock phase 

/ \ i P ~> 1<»<V2 ^ 

^ = 1 r /o , i <- <- r (3) 

I p + , L/2 + 1 < n < L. 

We then use the conservation of current to compare the current at two different sites, 
one at the impurity and the other far in the bulk of the low- or high-density domain. 
In the short-range model this gives densities p_ = r/(l + r) and p + = 1/(1 + r), where 
< r < 1 is the reduced hopping rate characterizing the defect there. In the long- 
range case, the conserved quantity is the site-independent particle current defined as a 
summation over all the particles jumping over and through some site k 

k L L L 

J = Y1 Yl Pi( r n{l - r n+l )) + Yl Pi(r n {l-r n+l )). (4) 

n=l l=k+l-n n=k+2l=L-n+k+2 

In the bulk of the low-density domain, i.e. far away from the impurity (k = L/4), the 
current j is roughly Al(<t)p_(1 — p_). On the other hand, the current over the impurity 
(obtained by setting k = L) is given approximately by (Al(ct) — l)p + (l — p_). [In both 
cases we have neglected terms that arise from the jumps longer than L/2, but they 
vanish in the limit L — > oo anyway] Equating these two expressions gives p_ and p + 
for a very large system in terms of A = Ylh^i lPi 

Equation (JSJ) may be reduced to the form of the short-range model, p_ = r' / (r' + 1) and 
p + = l/(r' + 1), where the reduced hopping rate r' is given by r' — (A — 1)/A. One way 
to interpret this expression is as the ratio of the average hopping length over the defect 
site (effectively reduced by one due to the defect site), — l)pz, versus the average 
hopping length in the absence (or far from) defect, Yli^Pi- 

Compared to the results of Monte Carlo simulations, we find that the above result 
becomes accurate only for large values of a, a > 3 (figure EJ). Therefore, it is not 
surprising that equating p_(cr c ) = p+(cr c ) in the above expression to obtain a c gives 
trivial value <rf MF = 1 ruling out the shock-free phase completely. Such a crude 
approximation that does not allow spatial variation of density cannot produce non- 
trivial transition to the shock-free phase in the short-range TASEP as well, where it 
gives r s c MF = 1 [13]. 

In order to explain the nontrivial value of o~ c > 1 and to study the scaling properties 
of the profile, we have to consider the site-dependent version of the mean-field approach. 
Let us first renumber the lattice sites so that n = —K, . . . , K, L = IK with an impurity 
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Figure 5. A comparison of bulk densities p_ in the shock phase obtained by the simple 
mean-field approach (horizontal lines) and by Monte Carlo simulations [L — 6400, 
p = 1/2 and t = 10 8 MCS/site). 



at the site n = 0. Starting from the master equation, lattice equations for the time- 
dependent average density profile (r n )(t) read 

d L L 

-u( T n) = ^Pi(r n _,(l - T n )) - Pi( T n(l - r ™+;)), (6) 

1=1 1=1 

l^n l^L-n+l 

where periodic boundary conditions have been assumed, Tx +n = T-K+n-i- m the mean- 
field approximation where (T n T m ) — > (r n )(r m ), the lattice equations (jUJ) in the stationary 
limit d(r n )/dt — > reduce to a system of L nonlinear equations in L variables restricted 
to the condition Yln=i( T n) = N, which can be solved numerically. For this purpose, we 
used the HYBRD algorithm taken from the MINPAC library [|| for various a and system 
sizes L. In the shock phase, the results reproduce those of the Monte Carlo simulations, 
but one generally finds a narrower shock and a shift in the density profile in each of the 
domains (figure [H^)- The first result is not surprising, as the fluctuations in the position 
of the microscopic shock are usually neglected in the mean-field approximation. Indeed, 
the width of the shock obtained from the mean-field data seems to obey the L 1-cr / 2 
scaling compared to the L 1 / 3 of the Monte Carlo data. On the other hand, we have no 
explanation for the shift at this moment, but we note that it is observed in the short- 
range TASEP with a defect site as well [T3]. Nevertheless, we find the site-dependent 
corrections in the region dominated by the long-range correlations described correctly in 
the mean-field approximation (figure [6b). This becomes more apparent in the shock-free 
phase, where p_ and p + become 1/2 and the profiles coincide with the results of the 
Monte Carlo simulations (figure [7]). 

In consequence, the same scaling of the profile decay from the impurity observed 
in MC simulations is reproduced in the mean-field approximation as well. Similar 
behaviour was observed in the pure case of this model with open boundary conditions 

\ http : / / www . netlib .org/ minpack / 
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Figure 6. (a) A comparison of density profiles obtained in the mean-field 
approximation and by Monte Carlo simulations for a = 1.8 (L = 6400, p = 1/2 
and t = 10 8 MCS/site). (b) The site-dependent corrections to the bulk density (t £ / 4 ) 
obtained in the mean-field approximation coincide with the results of the Monte Carlo 
simulations. 
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Figure 7. (a) A comparison of density profiles obtained in the mean-field 
approximation and by Monte Carlo simulations for a = 1.1 (L = 6400, p = 1/2 
and t = 10 8 MCS/site). 



[22], and indicates that the long-range hopping has suppressed the importance of 
fluctuations to the point that the mean-field approximation produces the correct scaling 
exponent. This motivates us to analyze the leading power-law contributions in the limit 
of large L, and obtain the analytic expression for the exponent v. To this end we define 
the site-dependent corrections <ft n to the (unknown) bulk densities p_ and p + in the 
stationary regime (d(r n )/dt — > 0) 

, . (p+ + 0n, -K<n<0 

\T n ) = < (7) 

{ P- + <Pn, < n < K. 
For < n < K, the r.h.s. of ([6]) gives the following terms in powers of <p n 
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where the following notation has been used, A+ n = <p n+ i — <p n and Aj"^ n = </>„ — 
Although these equations cannot be solved analytically, one can estimate the relevance 
of each term by taking the lattice spacing a = 1/L as a small parameter and comparing 
the order of each term in the limit a ^ 0, na — > x < oo. First term in cjP n is of the 
order 0(a CT+1 ), because p n ~ L"^ +1 ' x~^ a+1 ' . For < n < L/A, the second term in 0° is 
positive and of the order 0(a <J ). To estimate the cf) l n terms, we assume the scaling form 
(p n = L~ v f-(n/ ' L) with an unknown exponent v and /_(n/L) < for < n < L/4. 
The second and the third terms in <fi\ are then of the order 0(a u+a+1 ) and — 0(a u+a ), 
respectively. Neglecting the higher-order (nonlinear) terms 4> n 4>m, the only nontrivial 
terms to estimate are 



J2 PAt<t>n~J2 PlA l 
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(9) 



where we introduced Ap_ = 1/2 — p_ to distinguish two cases, p_ = 1/2 (Ap_ = 0) and 
p_ 7^ 1/2 (Ap_ ^ 0). In the pure long-range hopping on the infinite lattice with density 
of particles p, these terms can be estimated in the continuous limit a — > 0, x = na, 



<p n -> <f>(x) 
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Ap 
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(10) 



where L> CT = -T(-o-)cos(-na /2) / ((a + 1), £> 2 = CO - l)/(2CO + l )) an d A a denotes 
fractional Laplacian (also referred to as the Riesz fractional derivative [21]) defined 
with respect to the Fourier transform jF{A CT /(a;)} = — \k\ a f(k), f{k) = T{f{x)}. For 
1 < a < 2, the leading term in ([TO]) refers to the space-fractional diffusion A a <p, typically 
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found in the anomalous diffusion mediated by Levy flights (for a recent review see [25] ) • 
On the other hand, the leading contribution to the drift term (jlip is local for all o > 1, 
but the a-dependence is present in the collective velocity (1 — 2p)\(o~). 

Although the presence of impurity breaks the translational invariance and the 
continuous limit no longer holds, we expect that the estimation of the leading terms in 
f fTOj) and pip is still correct and given by 0(a a+u ) and -0(a u+1 ), respectively. Finally, 
turning back to the r.h.s. of ([7]) and taking the two lowest-order terms with opposite 
sign, the second term in 0° and the first term in ( fill) , we obtain the exponent v = a — 1 
for all a for which p_ ^ p + (i.e for a > er c ). For a c < a < 2, this result is exactly the 
same as the one conjectured from the Monte Carlo simulations, but it breaks down for 
a > 2 where the correlations become important and v takes the constant value v — 1. 

The shock- free phase is obtained by inserting p_ = p + = 1/2 in (jHJ), which gives 
only four terms 



\(Pr, 



PL-n+l) 
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+ 4>l(Pn - PL-n+1 ) +<Pn 
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(12) 



Since the drift term vanishes, the nonlinear effects become important and the density 
profile may have more complex form than the one given by the power-law behaviour 
with a single exponent v. This is indeed observed in the results of the Monte Carlo 
simulations, for which the assumed scaling relation <fr n = L~ u f(n/L) no longer holds. 

So far we established that the above mean-field analysis predicts the leading 
contribution to the exponent u(a) as long as a c < o < 2, but it gives no analytical 
expression for the bulk densities p_ and p + which would enable us to calculate the 
transition point o~ c . To this end we use a different argument which relies upon the 
scaling properties. If we suppose that there is no discontinuity in the exponent v(cr) 
at the transition point a c and use the result, argued below, that v = 1/3 at the onset 
of the shock-free regime, a c is obtained by equating u(a c ) = o c — 1 = 1/3. This gives 
a c = 4/3, which is also within the bounds 1.32 < a c < 1.4 estimated from the Monte 
Carlo simulations on large system sizes. 

Basically, the argument for v(o~ c ) = 1/3 is the same as the one that gives v — 1/3 in 
the shock-free phase of the short-range model |18j. First, one notices that the processing 
of density fluctuations across the impurity is reduced only for positive fluctuations, 
leading to the excess of particles in front and depletion of particles behind the impurity. 
If the centre of mass of density fluctuations travels faster than it spreads, one can 
ignore the spreading and estimate the number of excess particles solely by looking at 
the number of such fluctuations present in the system at some given moment. This 
number, however, cannot exceed the time tf needed for the fluctuation to travel across 
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the system. Since on average the positive and the negative fluctuations cancel out, the 

1/2 

number of excess particles 5N scales with tf as 5N ~ tj . Therefore, we can estimate 

the number of excess particles by estimating tf. In the shock-free phase, the centre- 

of-mass velocity of density fluctuations v g = Sj/5p ~ A[l — 2p(x)] is zero in the bulk, 

but remains nonzero near the impurity, v g ~ x~ v . The centre-of-mass position thus 

depends on t as xqm ~ t~ and so the time tf for the fluctuation to travel across 

the system scales with L as t/ ~ L u+l . The excess of particles then scales with L as 

5N ~ t l j 2 ~ L^ +1 ^ 2 and the self-consistency implies that 5N ~ J dxx~ v ~ L x ~ v , which 

gives u(a c ) = 1/3. This result is valid only if the spreading of density fluctuations 

1 / i 
£ ~ t /z is less than xcm ~ i.e. if v < z — 1, where z is the dynamical exponent. In 

the long-range case, the exponent z takes the value z = min{a, 3/2} [26], which means 

that cr c is also the limiting value at which the spreading becomes relevant. [For a > o c 

the spreading is irrelevant since the collective velocity v g is constant and nonzero in the 

bulk.] 

4. Conclusion 

In this work we study the asymmetric exclusion process with long-range hopping on 
the periodic lattice with single impurity site inaccessible to other particles and acting 
as a defect site. The key advantage of such impurity is that it obstructs the flow 
of particles but does not introduce any new parameter besides a. In spite of the 
long range of hopping, we show that for a > a c the impurity induces a macroscopic 
shock in the stationary state, similarly as in the standard TASEP model with a defect 
site. In particular, our model reproduces two main features of the short-range TASEP 
with defect site: namely, (a) the shock width scales with the system size L as 

£1/2 

or 

L 1 / 3 depending on the value of density, p^l/2orp = l/2 respectively, and (b) the 
density profile in the shock phase (a > o~ c ) displays a power-law decay away from the 
impurity. As in the pure long-range case without a defect, we show that the latter 
is well reproduced by the mean-field approximation and within this approximation we 
obtain the analytical expression for the a-dependent exponent of the power-law decay 
v — a — 1. 

Finally, by extending to the present model the argument of [18] that gives v = 1/3 
in the shock-free phase of the standard TASEP with a defect site and assuming the 
continuity of the exponent v{o~) at the transition point o~ c , we conjecture on the possible 
exact value of o~ c = 4/3 and find it to be in a good agreement with our numerical 
simulations. 

Although the model is not in the one-to-one correspondence with the standard 
TASEP with a defect site but rather shares some of its features, we believe that a 
clear example of the phase transition taking place at the nontrivial value of the control 
parameter, as in the presented case, may be relevant to the recent discussions on 
whether the static blockage in one-dimensional driven diffusive systems always induces 
a macroscopic shock [181 EZ]- 
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